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Abstract 

Sliding motion is evolution on a switching manifold of a discontinuous, piecewise-smooth 
system of ordinary differential equations. In this paper we quantitatively study the effects 
of small-amplitude, additive, white Gaussian noise on stable sliding motion. For equations 
that are static in directions parallel to the switching manifold, the distance of orbits from the 
switching manifold approaches a quasi-steady-state density. From this density we calculate 
the mean and variance for the near sliding solution. Numerical results of a relay control 
system reveal that the noise may significantly affect the period and amplitude of periodic 
solutions with sliding segments. 

1 Introduction 

Nonsmoothness and noise are two features of a dynamical system that may be the cause of 
important qualitative behaviour. Hybrid and piecewise-smooth systems are utilized in a wide 
variety of fields to model phenomena that involve switching, impacts or other nonsmooth elements 
[H [21 131 m ISl |6] . Recent studies have successively explained novel behaviour that may occur in 
such systems. For instance, so-called discontinuity maps have been developed and analyzed in 
order to understand border-collision scenarios in piecewise-smooth systems of ordinary differential 
equations El E]. Moreover, parameter uncertainty, background vibrations and other sources 
of noise are ubiquitous in real systems. Studies of stochastic differential equations have, for 
example, led to an understanding of noise-induced dynamics such as stochastic resonance and 
coherence resonance in excitable systems [TOl [HI [12l [13] . However, investigations into systems 
that are both piecewise-smooth and involve noise are relatively uncommon. 

One-dimensional, piecewise-linear maps with noise have been the subject of some isolated 
investigations [HI [151 [SI [13 [IB]- In [S], noise- induced transitions for a two-dimensional. 
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piecewise-smooth system of ordinary differential equations are explained through an analysis 
of a one-dimensional, piecewise-linear return map with additive noise. Simple vibro-impacting 
systems have been analyzed with stochastic averaging ^U\. Exact results are attainable for a 
classical, unforced, linear oscillator with elastic impacts [211 [22], whereas more complex scenarios 
have been investigated asymptotically and numerically [231 [21]. Stochasticity in switched control 
systems is particularly important in regards to robustness of an output signal to noise [23 [211 [2Z] . 
Here Lyapunov functions are an invaluable mathematical tool for determining stability because 
they do not necessitate smoothness in the vector field. Noise-induced oscillatory motion has 
been studied in piecewise-linear systems for which local linearity makes some key calculations 
tractable [281129] . 

In this paper we consider piecewise-smooth, stochastic differential equations for which the 
underlying deterministic dynamics are described by the ODE 

X = Fi(x) , ^eQi , (1) 

where each Qi G V G is open, nonempty and pairwise-disjoint, Uifli = V, and each 
Fj : f2j — )■ is a smooth function. Equation ([T]) is a Filippov system [30] and well-suited 
to model phenomena that alternate between different dynamical regimes, such as vibrating sys- 
tems experiencing impacts or friction [311 [321 |33l [311 [35], and switching in electrical circuits 

[2 [3 EE]. 

Boundaries between the neighbouring subdo mains, Qi, are codimension-one surfaces termed 
switching manifolds. Often, a section of a switching manifold has the property that on either 
side of the manifold the vector field points towards the manifold. In this case any orbit that 
reaches the switching manifold becomes trapped on the manifold for some time. The resulting 
motion on the switching manifold is known as sliding motion, Fig. [Tj-A. Formally this is achieved 
by Filippov's method [371 [301 [I] which defines a vector field on the switching manifold by the 
unique convex combination of the two limiting vector fields on either side that is tangent to 
the switching manifold. Sliding motion corresponds to the sticking phase of stick-slip oscillators 
[381 [39] and the coalesced regime of a piecewise-linear relay control system [lOl [H] . The addition 
of small noise pushes orbits off the switching manifold, but large excursions are curbed by the 
deterministic component of the system. Fig. [T]-B. Thus the motion is balanced by the competing 
actions of noise and drift. In this paper we explore these dynamics more carefully. 

We here summarize the basic effect of adding small noise to a smooth system. Below we 
compare this to our results for the piecewise-smooth system ([T]). Consider an A^-dimensional 
system 

X = $(x) , (2) 

where $ is a smooth function. The addition of small amplitude, time- independent, white, Gaus- 
sian noise gives the stochastic differential equation 

d^ = <!>{^)dt + y/eB{x)dW{t) , (3) 

where < 5 ^ 1, B{x.) is an N x N matrix with a smooth dependency on x, and W{t) is a 
standard Brownian motion [121 [13]. Let ^^(x, t|xo) denote the transitional probability density 
function (PDF) for the point x(t), given x(0) = Xq. A straight-forward expansion in powers of 
^/e reveals that the mean of ^^(x, t|xo) differs from the deterministic solution (the solution to ([2]) 
with x(0) = Xq) by 0{e), whereas deviations are 0{^/e) [131 [11]. 
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Figure 1: Schematics of a Filippov system near a switching manifold that attracts orbits from 
both sides in the absence of noise, panel A, and with small amplitude additive noise, panel B. 



In this paper we derive an analogous result for sliding motion for which the above method 
of expansion does not work because the vector field is discontinuous. Instead we analyze a 
one- dimensional, discontinuous stochastic differential equation for a quantity representing the 
distance from the switching manifold. We find that the mean solution differs from Filippov's 
deterministic sliding solution by 0{e), and deviations are 0{y/e), matching the smooth case. 
Moreover, the calculations suggest conditions necessary for the noise to induce a change in the 
dynamics that is not dominated by randomness. 

Here we outline the remainder of the paper. We motivate our work in ^ by illustrating 
that noise may significantly reduce the amplitude and period of a solution to a prototypical 
relay control model involving segments of sliding motion. In §3] we introduce a system of two- 
dimensional stochastic differential equations, fl5])- ffTT]) . that describes stochastically perturbed 
sliding motion relating to a linear switching manifold in the case that system is the same in 
directions tangent to the switching manifold. In this case, the equation for motion in the direction 
orthogonal to the switching manifold, x, is independent of the variable representing displacement 
tangent to the switching manifold, y, and for this reason is amenable to an exact analysis. We 
leave a description of more general scenarios for subsequent work. In [|l]we analyze the stochastic 
differential equation for x{t) and derive its quasi-steady-state distribution. In [|5]we analyze the 
equation for y{t) and obtain expressions for the mean and variance of y{t). Derivations for this 
section are given in ^ and Appendix [Al Conclusions are presented in ^ 

2 Periodic orbits with sliding and relay control 

Periodic orbits involving sliding have recently been described in models of relay control systems. 
Broadly speaking, a relay control system is a system that aims to control a variable using the 
measurements of an input signal via a switching action HHl SZl HE]- Relay control systems 
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are commonly modelled by 



X = y4x + Bu , 

^ = C"^x , (4) 

u = -sgn(v9) , 

where x G M^, (p is the signal measurement and u is the control response, [101 El [1]. The system 
(jl]) is a Filippov system with a single switching manifold, {x | C"'"x = 0}, on which sliding 
may occur. In this system sliding corresponds to the idealized scenario of discrete switching 
events occurring continuously in time. Periodic orbits of (jl]) that involve sliding are described in 

(nmsiEoiEi]. 

As an example we consider the following canonical form, taken from [H] (also given in [1]), 
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where C ^ is a parameter. Fig. [2] illustrates a stable periodic orbit of @ with ([5]) involving 
sliding motion on the switching manifold. There are 12 separate sliding segments per period. 
These correspond to time intervals for which xi (the first component of the vector, x) is zero. The 
stability of periodic orbits with sliding may be determined by analyzing the Jacobian of a return 
map [STl [52] . The robustness of periodic orbits with sliding has been briefly investigated by 
studying the size of the basin of attraction of the periodic orbit [52], and imposing a short time 
between consecutive switching events [53] . 

Fig. [3]-A shows a typical orbit for the system when noise is added to the control signal as 

dx = {Ax-B sgn (C"^x)) dt + ^/^B dW{t) . (6) 

By comparing with Fig. |2]-B, this shows that the noise may dampen the oscillations and induce 
an increase in frequency. Fig. [3]-B shows that the average length of time per oscillation of 
decreases as the noise amplitude increases. We have observed similar behaviour for different 
values of (. In §7] we use the results below to speculate on the reason for this behaviour. 



3 A simple set of equations for stochastically perturbed 
sliding motion 

We are interested in the effects of noise on the dynamics of ([1]) near a switching manifold. In this 
section we introduce a simple set of equations that approximates ([1]) near a switching manifold, 
then formulate the inclusion of noise. 

If a switching manifold of ([T]) is locally we can choose a coordinate system such that in 
a neighbourhood of the origin, x = 0, the switching manifold is simply e^x = 0{K) [9]. In this 
paper we are not concerned with effects due to nonsmoothness in the switching manifold and for 
this reason suppose that the switching manifold is the coordinate plane, e^x = 0. Two smooth 
subsystems govern the nearby flow, thus, locally, we may write the deterministic system as 

r FW(x), e7x<0 . 
^ - I F(R)(x) , e7x>0 ' 
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Figure 2: A stable periodic orbit of dl]) with (E]) and C = —0.06. Here x = (xi, a;2, Xs)""". The 
periodic orbit exhibits shding on the switching manifold, Xi = 0. 



where F^^^ and F^^'^ are, say, C^. 

Suppose there exists a section of the switching manifold, call it S, for which e7-^*-'^''(x) > 
and e7F(^)(x) < 0, as in Fig. [T]-A. In this scenario, forward orbits arrive at S from either side. 
We use Filippov's definition to define dynamics constrained to S ^S7\ [501 [T]. S is known as an 
attracting sliding region and evolution on S is referred to as sliding motion. 

However, with additive noise, the system ([7]) seems to be too complex for us analyze to a 
degree of detail that is useful. This is because both F*^^) and F^^"* depend on all components of 
the vector x, and we have been unable to analytically solve the resulting A^- dimensional stochastic 
differential equation. Consequently, for this paper, which represents a first detailed analysis of 
sliding motion with noise, we ignore the dependency of F^^^ and F*^^-* on components of x parallel 
to S as this enables us to reduce mathematical problems to one dimension but still capture what 
seems to be the essence of stochastically perturbed sliding motion. Moreover, this provides a 
useful approximation to the general case over short time-frames. 
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Figure 3: Panel A shows a time series of (j6]) with e = 0.001, using the same parameter values as 
in Fig. |2j Panel B shows the median, upper quartile and lower quartile values of 1000 numerically 
computed oscillation times, for several values of e. To obtain each oscillation time, we computed 
an orbit up to t = 100 and identified the last three instances at which the value of changed sign 
(discounting rapid sign changes over a handful of grid points near this value due to noise), then 
subtracted the first time from the third time. Orbits were computed with the Euler-Maruyama 
method of fixed step size. At = 0.0001. 



Given that F^^^ and F^^^ are functions of only e^x, the remaining — 1 components of x 
may be treated identically and for this reason it suffices to study a two-dimensional system. We 
let X = [a; y]""" and add small amplitude, white, Gaussian noise independent to the state of the 
system. Assuming for simplicity that the noise in x is independent of the noise in y, the resulting 
stochastic differential equation may be written as 



dx 
dy 



(f){x) 
ip{x) 



dt + y/e 



dWi{t) 

KdW2it) 



(8) 



where Wi{t) and W2it) are independent Brownian motions, < e ^ 1 and k > are constants, 



and (f) and if) are piecewise-C^ that for small |x| are given by 



j.( \ _ { ai + clx + o{\x\) , x<0 , . 

'^^''^ \ -aR + CRX + o{\x\) , x>0 ' ^ ' 

Hx) = + + • (10) 

^ \ Br + dRX + o{\x\) , x>0 ^ ^ 



We assume 



aL,aR>0, (11) 

to ensure that in the absence of noise the switching manifold {x = 0) is an attracting sliding 
region. Since (j) and ifj are independent of e, their coefficients are 0(1). Consequently, for x(0) 
near zero, orbits of ([S]) likely remain near x = for relatively long periods of time, as shown in 
§4.H and for this reason we do not specify the behaviour of and ip for large \x\. 



4 Properties of x(t) 

Since (E]) lacks dependency on y, the equation for dx is decoupled from y: 

dx = (j){x) dt + ^/6 dWi{t) 

Given x(0) = xo, let Pe{x,t\xo) denote the transitional PDF for the value of x{t), as governed by 
( fT2l) . Despite the discontinuity at a; = 0, Pe{x, t\xo) is unique and continuous on M x (0, oo) x M. 



For x ^ and t > 0, the PDF satisfies the Fokker-Planck equation 

dt dx 2 dx^ ' ^ ^ 

with the initial condition p^{x,0\xq) = 6{x — Xq), [l2l HSl |5l]. In ||6] we provide an explicit 
expression for pe{x,t\xo) in the special case that (p is piecewise-constant. For the remainder of 
this section we use (fT3l) to determine the long time behaviour of x{t). 

If 4>{x) > for all X < 0, and <f>{x) < for all x > 0, then ( |T2|) has a steady-state density on 
M centered about the origin. Otherwise, (p{x) = for some x ^ 0, and with nonzero probability 
orbits may cross this value of x and undergo dynamics far from the origin not described by the 
expansion (Q. However, regardless of the global nature of (p, since e is small the local attraction 
to the origin is relatively strong. Thus we expect orbits to remain near the origin for long periods 
of time and be distributed by a quasi-steady-state distribution for large but finite t. In §4.11 we 
determine the mean escape time of orbits from an 0(1) neighbourhood of the origin. In §4.21 we 
use f|T3l) to derive the quasi-steady-state probabihty density function asymptotically. 

4.1 Escape from a neighbourhood of x = 

For the function (j){x), iQ, with ffTTl) . in the case that (p{x) = for some a: 7^ 0, it is necessary to 
identify a value Xb > 0, independent of e, such that 

min \<j){x)\ > ^mm{aL,aR) . (14) 

|a;|<i'i, Z 
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Then in the interval [—Xb,Xb], the drift of f|T2|) is towards x = 0. With e > and any initial 
condition xq G {—Xb,Xb), x(t) will eventually escape [—Xb,Xb] with probability 1. Calculating the 
time to escape is a standard problem in the context of a single potential well, where the potential 
function is given by 

U{x) = - [ <P{y)dy. (15) 



The mean escape time, T{xo), may be found exactly SH US] . Via Laplace's method of asymp- 
totic evaluation of integrals |56j, it follows that whenever |xo| < Xb there exist ^-independent 
constants ai and a2 such that 

— — 

T{xq) ~ eaie~ . (16) 
For instance if U{xb) < U{—Xb) (other cases are similar), 

T(.o) ~ ^(^£±^e^ . (17) 
^ °^ 2aLaRU'ixb) ^ ' 

4.2 The quasi- steady- state probability density function 

In §4.11 we showed that the mean escape time from an e-independent neighbourhood [— Xf,,x;,] 
is exponentially large in K Consequently, we can assume that the probability an orbit escapes 
[— within the polynomial time -pj, for any fixed M > 0, is extremely small. We let 

i = e^t , (18) 

represent the long time scale, and look for a solution to the Fokker-Planck equation f[T51) as a 
function of x and t. By substituting (fT8l) with the WKB-type expansion |l2l H3] . 

Pe{x,t\xQ) = e J , (19) 



into ffT3l). we arrive at 



j^dqe _ ( dqe d \ [1 dqe 



Therefore 



dx 

and so by integrating 0(x) we obtain 

= Tit) + 1 ! t "'^.k ! ? ° K o(^«) . (22) 



—2ajix + crx + o(x ) , X > 

for an e-independent function r. As a function of a; and t, the dependence of this solution on t 
and appears only in 0{e^) terms which may be ignored. Consequently we treat the solution 
as solely a function of x and refer to it as the quasi- steady- state solution, Pqss,e{.x)- Specifically 
(fT9|) and (!22|) combine to give 

r :^g^(2aiX+cz,x2+o(x2))+0{e*^) x < 

Pqss,e(a;) = < £.^i(_2a^x+CBx2+o(x2))+O(£A'0 3; > ' ^^^^ 



where we must have 

K,_^^-J^£!i±p^. + o(e), (24) 

ai + ttR aittRlaL + aRy 

to ensure Pqss,e is normahzed. 

For small e and xq, the transitional PDF of f|T2l) . p£(x,t|xo), quickly settles to fl23l) . The 
scaling 

X ^ t , ^ 

a; = - , ^ = ~ ' (25) 

transforms (fT2|) to 

rfx = rft + diyi(t) , (26) 

from which we infer that Pe{x,t\xQ) approaches (l23l) on an 0{e) time-scale, when xq = 0{e). 
Furthermore, for times in the range e^~^ < t < £~*^, where 5 > 0, it is reasonable to suppose 
X ~ Pqss,e, in which case 

{xsgn{x)) = + . ■; + 0(.'). (29) 

2aLaR{aL + aR) 

which are useful in the next section. 



5 Moments of y{t) 

In this section we compute the mean of y{t) and conjecture the leading order term of its variance. 
We assume x ~ Pqss,e at all times under consideration which greatly simplifies calculations. We 
begin by deriving y{t) when e = 0. 



5.1 Deterministic sliding motion 

When e = 0, ([8]) is the Filippov system: 



X 










i/j^x) 



(30) 



As in [271 [201 [I], we define a vector field for sliding motion on the switching manifold {x = 0) 
by the unique convex combination of the two vector fields at the manifold that is tangent to the 
manifold. That is, 

X 

y 



q) 



slide 



hL 



+ q 



-aR 
bR 



for the unique scalar, q G (0,1), for which Xsiide 
therefore 

_ aRbL + aLbR 

?/slide — ; 

ai + aR 



0. Solving Xsiide = gives q 



(31) 
and 



(32) 
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Consequently, if (x(0),y(0)) = {0,yo), then Xsiide(i) = and 

yshde{t) =yo^ : 1 ■ (33) 

5.2 The mean of y{t) 
From dH]) and (11 01) we have 

dy = {^-^-''-^s^n{x{t)) + '^^ 

+ v^dM^2(t) . (34) 
Integration yields 

yif) = yo+ ^^^ ~ 2 / sgn(x(g)) + 2 / ^^^^ ^'^ 



d — d /"* 

/ a;(s) sgn(x(s)) c?s + / o{\x{s)\) ds + ^/eKW2{t) 

2 Jo 



and therefore 







t j-t 
(x(s) sgn(x(s))) (is + / (^o{\x{s)\)') ds . 



{y(t)) = y« + '^t-''J^j\.!,.(M))ds + 'i^j\.(s))ds (36) 

2 

If a; ~ Pqss,e, by substituting ( I27|) - fl29|) we obtain 

(vit)) = v....{t) + 1"^'''' - " + °'^> ■ 

where the ^-independent terms have combined to form ysWAeif)^ fl33l) . 

Therefore as e — )■ 0, the mean of y{t) hmits on Filippov's sliding solution, ysVidcif)- This is 
non-trivial because Filippov's method, to obtain (l33l) . and standard stochastic dynamical systems 
definitions, to obtain (137|) . are not immediately related. Note that the perturbation of {y{t)) 
from ysiide(^) is order £, mirroring the result for smooth systems, see The explicit expression 
for the coefficient of the 0(e)-term in (l37Il is particularly useful. For instance, we can see that if 
cl = cr and di = dR, then we require the asymmetry gr in order for the 0(e)-term to be 
nonzero. 

5.3 The variance of y{t) 

The variance of y{t) may be computed via 

Var(y(t)) = {y{tr) - {y{t))' , (38) 
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however this requires knowledge of J9e(x, t\xQ), for which we have not been able to obtain a useful 
expression in the case of general (p. We conjecture that the leading order terms of Var(y(t)) are 
independent of non-constant terms in (p and ip because x{t) = 0{e) with high probability. Indeed 
this is consistent with numerical simulations, Fig. HJ In view of the result for the case that (j) and 
ip are piecewise-constant (Theorem [21 given below), we propose the following result: 

Conjecture 1. Consider ^ with / flT]) . Suppose x{0) is random with PDF, Pqss,e, andy{0) = y^. 
Then for any 5 > 0, whenever e^~^ <t< C^'^ we have 



{hL - b 



\2 



Var{y{t))=eKt + ^ ^et + Oie'). (39) 

From here until the concluding section, ^ we study the case that (j) and are piecewise- 
constant, i.e. 

= (^"' "^n ^ (40) 

' \ -aR , x > ' ^ ^ 

, / N ( br , X < , , ^ . 



Theorem 2. Consider ^ with 07]) and suppose (p and ip are given by ^^D^-^^^. Suppose x(0) 



is random with PDF, Pqss,£, and y{0) = yo. Then for any 5 > 0, whenever t > e we have 



{hL - b 



Var{y{t))=eKt + P^ — "-^^et + 0{e') . (42) 



\2 



(ttL + aR) 

We prove this result in the next section. 



6 Two- valued drift and a proof of Theorem [2 



The stochastic differential equation f[T2|) with fl40|) : 



dx 



-aR 



X <0 
X > 



dt + y/edWi{t) , 



(43) 



Pe{x,t\Xo 



has been referred to as Brownian motion with two- valued drift. The transitional PDF of this 
process was first derived by Karatzas and Shreve in [57]. In this section we state this PDF and 
use it prove Theorem O 



6.1 The transitional probability density function for x{t) 

The transitional PDF for ( l43l) is given by 

h^{t,b,aR) *h^{t,b-x- Xo,aL) + G'absorb,e(a:, t, ailxo) , 
" hs{t, b + x, Or) * hs{t, b - xo, ai) db , 



e „ 

2 —2af}X 



he(t, b + Xq, aR) * h^it, b-x, ai) db , 

e^^ Jo hs{t,b + X + Xq, Gr) * he(t, b, ai) db + G'absorb,£(2;, t, -aR\xQ) 



Xo < 0, X < 
xo < 0, X > 

Xq > 0, X < 

xo > 0, X > 
(44) 
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Figure 4: The six horizontal bars are 95% confidence intervals for Var(y(t)) — + (a^_^a^)2 j 
when (a^, a^, bi, &_r) = (2, 1, 1, 0), k = 0, £ = 0.01, t = 1, and and -0 are piecewise-linear using 
different values of cl, cr, di and dR as indicated by the table. Each confidence interval was 
obtained from 10® orbits computed numerically with the Euler-Maruyama method of fixed step 
size, At = 0.0001. The results are consistent with Conjectured] which predicts that the values of 
Cl, cr, di and d^ affect the magnitude of Var(y(t)) by at most 



where 



, , , \Xo\ i^o-f^tf 

he{t,xo,fi) = -j===e , (45) 



is the PDF for the first passage time to zero of Brownian motion with constant drift, 

<Jabsorb,£(a:,t,^|xo) = ^ e -e - e 2.* . (46) 

V27r£t \/27ret 

is the transitional PDF for Brownian motion with constant drift and an absorbing boundary 
condition at zero, and 

m*Mt)^ f h{r)h{t-r)dT , (47) 

is the convolution relating to Laplace transforms. Fig. [5] shows (jSj) at different times. 

In [57] , Karatzas and Shreve derive (jH]) for e = 1 by using Girsanov's theorem [58] [59], [60] 
and the trivariate PDF of Brownian motion, its positive occupation time, and its local time about 
zero. The result for e ^ 1 follows simply from the scaling fj25]) . The Laplace transform of (jHJ) 
can be written as an integral- free expression and this was achieved in the earlier paper [61]. We 
do this below for xq = and x > 0. In [521 IHS], the authors derive Pe(0,t|xo) and use this 
to bound PDFs for a large class of scalar stochastic differential equations. In [M], Pe{x,t\Q) is 
studied in the case aL,aR < 0. In [65j, Zhang derived an expression for the transitional PDF of 
Brownian motion with a general bounded piecewise-continuous drift function. This could be used 
to analyze the PDF of (fT2l) with general (j) asymptotically, but such a calculation is beyond of 
scope of this paper. 
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Figure 5: The probability density function of x{t), (jSj), with xq = 0.02, aL = 2, an = 1 and 
e = 0.01 at four different times. 



6.2 Proof of Theorem [2] 

In the case of two- valued drift (HOll . the quasi-steady-state density (!23ll is a true steady-state 
defined for all a; G M: 

Pss,e X = < ^ K = — . (48) 

[ Y^^^ , X >0 ttL + ttR 

With the notation 

we have the following expression for the probabihty that x{t) > 0, given xq = 0. 
Lemma 3. For any t > 0, 

p,{x,t\0)dx = -- J^aRPeiO,s\0) + --^{0,s\0)ds. (50) 



Proof. 

For X > 0, in the case of two- valued drift (l40l) . the Fokker-Planck equation (fT3|) is 



dpe d [ edps , 



Integration over x G [A, oo) yields 



p,(x,t|0)c?x = -(^a«p,(A,t|0) + |^(A,s|0)) , (52) 



for any A > 0. Then integrating with respect to t and taking A — )■ produces (!50|) where we 
also use 

lim / p^{x,t\0) dx = - , (53) 
i^o+ Jo 2 
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which may be demonstrated by noting that as t — 0+, Pe{x,t\0) is well-approximated by a 
zero-mean Gaussian. □ 
Proofs of the following two lemmas are given in Appendix |X1 Theorem |2] is an immediate 
consequence of Lemma M combined with (^7^ and 



Lemma 4. For the density 

f a*(0,*|0) + iM(0.<|0)* . (54) 

t [a^MO, m) + |||(0. tm) dt ^ ^^^^^^^ . (55) 



Lemma 5. Consider and suppose ai, > and x{0) is random with PDF, Pss,e- Then for 
any S > 0, if t > e^^^ , we have 

[ [ {sgn{x{s)) sgn{x{u))) ds du = l +- — '^^^ + 0{e^) . (56) 

Jo Jo [a-L + (IrY [ai + dRr 



In the special case, ai = = a, we can write Var{y{t)) exactly. In this case, by symmetry, 
Pi;{x,t\0) dx = |, thus by Lemma [3l ape{0,t\0) + |^-(0,t|0) = 0. Consequently the integral 
on the right-hand side of ( 175]1 vanishes, and ( ISTj) leads to 



Var(y(t)) = 3--! + ^ -T-^-^ e 




a t 

2e 



-.-'4+t' + ^]erfc[^] 1 , (57) 



which is consistent with (HI 



7 Conclusions 

When small noise is added to a Filippov system, orbits no longer slide along an attracting sliding 
section of a switching manifold. Instead, with high probability, orbits follow a random path near 
the switching manifold. Fig. [1]-B. The average size of deviations from the switching manifold is 
governed by the strength of the noise relative to the magnitude of the vector field in a direction 
orthogonal to the switching manifold. 

The general A^-dimensional stochastic differential equation formed by adding noise to the 
Filippov system ([7]) is particularly difficult to analyze due to the discontinuity and multiple 
dimensions. For this reason we made the supposition that the system is invariant along the 
switching manifold. This prevents an exploration of the effects of noise on orbits that reach 
the end of an attracting sliding region, but enables calculations to be reduced to one dimen- 
sion. Moreover, with this reduction directions parallel to the switching manifold may be treated 
identically and hence it is sufficient to study the two-dimensional system ([8]). 
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For the system ([8]), x{t) denotes the displacement from the switching manifold and is governed 
by f lT2|) with ([9]). Sample paths of f|T2|) settle to the quasi-steady-state PDF, Pqss,e fl23|) . on an 
0{e) time scale. This PDF is not a true steady-state because orbits escape a neighbourhood of 
zero. However, since this occurs on an exponentially long time scale, see §4.1^ it is suitable to 
assume x{t) is distributed by Pqss,e at times in a given range e^~^ < t < e~^^ , for any 5, M > 0, 
where we may take small 6 and large M, such that this is a long time interval. This assumption 
has the benefit of significantly simplifying our calculations. 

Equation ( l34ll is the stochastic differential equation for y{t), which represents displacement 
along the switching manifold. In the limit, e — )■ 0, the mean of y{t), (|371) . limits on Filippov's 
sliding solution, ?/siide(0- £ 7^ 0, the perturbation of from ysudeit) is 0{e), as for a 

generic smooth system. The perturbation depends on linear terms in and i/j. In order to gauge 
the effect of this perturbation on overall dynamics, it is necessary to compare it to the standard 
deviation of y{t). In the case that (p and i/j are piecewise- const ant, Var(y(t)) is 0{e), see Theorem 
HJ For general and i/j, we conjectured that the leading order terms of Var(y(t)) are unchanged. 
Consequently deviations of y{t) from are 0{y/e). Therefore, assuming e: <^ 1, we expect 

these deviations to dominate the difference between {y{t)) and ysiide{t)- 

Although our above calculations are for differential equations that are static along the switch- 
ing manifold, we can apply the basic principles gained to more general systems such as the relay 
control system, (jl]) with First note that since the deterministic equations are independent 
of e, when we make statements asymptotically in e, we implicitly assume that the magnitude of 
any parameter in the deterministic equations is much less than, in particular, But, for (jl]) 

with (j5]), = — C — 100, which is an extremely large value (we used ( = 0.06). Since xi rep- 
resents displacement from the switching manifold and X2 is a direction parallel to the switching 
manifold, in the context of (l8l)-( |T0|) . ^ corresponds to the values of cIl and d^. These values 
influence, to lowest order, the perturbation of the mean from the deterministic solution, ( 137|) . 
but not the deviation of sample paths from the mean, (1391) . This is a possible explanation for the 
noise-induced effect identified in ^ Specifically we found that when the parameters of (jlj) with 

are tuned such that in the absence of noise orbits settle to oscillatory motion with sliding, 
the addition of noise may significantly decrease the average oscillation time. Fig. [21 A further 
analysis is required in order to make more definitive statements. We believe that the decrease 
in oscillation time with e is a consequence of particular geometrical orientations and that an 
increase in oscillation time with e is equally possible for this type of system. 

This paper leaves many avenues for future investigations. Perhaps foremost, we would like to 
understand the perturbation of {y{t)) from ysnde in the case that the system explicitly depends 
on y. In particular we would like to understand how noise effects dynamics near the end of an 
attracting sliding region. Also it remains to consider more general forms for the noise, such as 
coloured noise or noise that is correlated in x and y, and study the effects of noise on other 
scenarios such as sliding bifurcations 
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A Proofs 



A.l Proof of Lemma [4] 

By Lemma Eland (iHD, 



which verifies (l54l) . To obtain (155|) . we multiply ( 152|1 by t and take A — i- to obtain 

t-J^ Pe{x,t\0)dx = -tiaRPe{0,t\0) + --^{0,t\0)] . 
Integration by parts yields 

^ t(^aRPeiO,t\0) + ^^{0,t\0)^ dt = ^ p,{x,t\0)-p,{x,T\0)dxdt., 
for any T > 0. Taking T — )■ oo gives 

\ poo POO 

Ps{x,t\0) - Pss,e{x) dx dt . 

We evaluate the double integral using Laplace transforms. For a; > 0, 



l^l^anpMm + l^mO)^ dt = 1^ 



roo _2 

P,,,,(x,A) = £b,s,.(2:)]= / e-^'Ke^dt 

Jo 

The Laplace transform of h^, is 



for A > 0. Using fj44l) and expanding about A = 0: 
P,(x,A|0) = C[peix,t\0)] 

POO 

/ H^{X,b + x,aR)H^{X,b,aL) db 
Jo 



2 D a: 

2e-^ 



2e^ 



-a/j + A/a|j + 2£A — a/, + A/a| + 2eA 



2e " aiaR ^ g(a| + a|) - 2a\aR{aL + Qfl)^: _^ ^^^^ 

+ aij)A 2aLaR{aL + aij)^ 



By dHD and ([MD, 



A|0) - P.,(x, A) = ^-K + 4)- 24«„K + aH). ^^ ^ 

eaiaRlai + any 
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and thus from fl6Tl) 



t[aRPe{0,m + ^^{0,t\0) ] dt 



£(a| + a|) - 2a\aii{aL + aR)x -^^r^ 



SttLaRlaL + an) 

-e{aL - an) 



dx 



2aLaR{aL + aR) 



as required. 

A. 2 Proof of Lemma [5] 

Since (fT2|) has no exphcit time-dependence, 



(66) 
□ 



t rt 



Jo 



(sgn(a;(s)) sgn(x(-u))) ds du 



/ (sgn(x(0)) sgn(x(|'u — (iw 
'0 Jo 

2 / (t — u)(sgn(x(0)) sgn(a;(tt))) (iu . 



(67) 



Furthermore, 



(sgn(x(0)) sgn( 



x[u 



CO roo 



sgn{x)pss,e{x) sgia{y)pe{y, u\x) dydx 



OO ^ — CXD 





Ke- 



oo 
oo 



+ / Ke- 
'o 



sgn{y)pe{y,u\x) dydx 

X) 

sgn{y)pe{y , u\x) dydx . 



(68) 



To evaluate (l67j) using fl68|) . we reorder the integrals of y, x and u, but to do this we must first 
transfer the w-dependence from the integrand to the limits of integration. We have 



sga{y)pe{y,u\x < 0) dy = 2/ p^{y , u\x < 0) dy - 1 



(ru roo \ 

J he{s,x,-aL) J pe{y,u - s\0) dydsj - 1 , (69) 



where in the second line we conditioned over the first passage time, s, of (|T2l) from x < to 0. 
By Lemma [3], and since K{s, x, —ai) ds = 1, for all a; < 0, we obtain 

/oo /•oo 
sgn{y)p^{y,u\x <0)dy = - hs{s,x, -ql) ds 
-oo J u 

-2 K{s,x,-aL) aRp,{0,v\0) + -^{0,v\0) dviSQ) 
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A similar calculation with a; > produces 

/oo poo 
sgn{y)ps{y,u\x > 0)dy = / h,{s,x,aR)ds 
-oo J u 

-2j h,{s,x,aR) J aRPe{0,v\0) + --^{0,v\0)dvd^71) 
By combining ( 1671) . (1681) . ( 1701) and (ITTl) we arrive at 

/ / (sgn(x(s)) sgn(x(M))) (is = 2 / {t — u) l Ke / hs{s, x, —ul) ds dx du 

Jo Jo Jo J-oo Ju 



fU — S 



+ 4 / {t — u) i Ke^~ / he{s,x,—aL) / aRPe{^-,v\Q) + -^^{^,v\{^) dv ds dx du 

Jo J~nn Jn Jo 2 (JX 



-oo 

o roo 



+ 2 / {t — u) I Ke e / he{s,x,aR) ds dx du 
Jo Jo Ju 

/■* /-^ -2atix r r-' edv^ 

-A {t-u) Ke^- / h,{s,x,aR) / aRPe{0,v\0) + -^{0,v\0) dv ds dx du . (72) 
Jo Jo Jo Jo 2 ox 

We now simplify the four terms of (!72|) . We define 

(^^(t, a)= / (t — u) / iCe e / hs{s, X, a) ds dx du , (73) 
Jo io io 

so that the first term of 0721) may be written as 

rt pO „ POO 

I I 2aj^x I CLr rj 

2 {t — u) Ke ^ / he{s,x, —ul) ds dx du = 1 —2Qe{t,aL). (74) 

Jo J-oo Ju «L + «R 

Similarly the third term is 

2 (t-u) Ke^^ dxdu-2Q,{t,aR) = — - 2Qe{t, ur) . (75) 

Jo Jo + 

To the second term of ( 1721) we reorder the integration such that dv is the outer-most integral 
instead of the inner- most integral. This step is straight-forward but requires some care with the 
limits of integration. We obtain 

4 / (t — u) / Ke s / hs{s,x,—aL) / aRp^{0,v\0) -\ —^{0,v\0)dvdsdxdu 

Jo J-oo Jo Jo 2 ox 

= 4^ (^aRp,{0,v\0) + ^^{0,v\0)^Q,{t-v,aL)dv , (76) 

and similarly for the last term of (|72l) : 

/•* /•°^ -2afl^ /•« /•"-'^ edp-^ 

— 4 / (t — u) / iTe £ / h^{s,x,aR) / aRPe(0,v|0)H — ^(0, t;|0) (if rfsfixc/M 

io io Jo Jo 2 ax 

= -4^ (^^p^{Q^v\Q)+^-^{S),v\Q)^Qe{t-v,aR)dv . (77) 



We are now able to write f l72|) as 

-t rt 



^0 



sgn(x(s)) sgn{x{u))) ds du = - 2{Qe(t, ul) + Qs{t,aR)) 

+ 4 ^ (^aRp,{0,v\0) + 1^(0, v\0)^ {Q,{t - v, a^) - Q,{t - v, a^)) dv . (78) 

Via OS]), it may be demonstrated that ajiPe{0,v\0) + |^^(0,f|0) decays exponentially to zero 
as f — >■ oo on an 0{e) time scale. For this reason it is helpful to apply the substitution v = - to 
obtain 



t pt 



^0 



sgn(x(s)) sgn{x{u))) ds du = - 2[Qs{t,aL) + Qe{t,aR)) 

+ 4£ (^ajtPeiO, ev\0) + |^(0, ev\0)^ {Qe(t - ev, a^) - Qeit - ev, an)) dv ,(79) 

and expand Qe{t — sv, ai) — Qe{t — ev, aji) in e such that the integral on the right-hand side of 
rgj) has the form 

Y.Y.^i,e' (ai?Pe(0,e5|0) + £M(o,e^|o)^ dv , (80) 



« 



for some coefficients aij. To obtain the coefficients, we first evaluate f l73p via multiple applications 
of integration by parts: 



3 , 5 



a t 



4a 4a3 Aa^ .pT^ I 6 3a2 2 



27r y 

, at^ et^ sH \ ^ ( a\ft 

K\ \ H erfc ' 



12 4a 4a3 4aV Vv^ 
In view of ( 17U]) . we use (ISTI) to obtain 

g.(t,«L) + g.(t,a^) = - TT^If^^ + 0(5^) , (82) 

2 2alaj^{aL + an) 

Q,(t-£^;,aL) -Q,(t-£v,a/j) = — — — + ^ ^ 2 / " — V + ^ — T ^ + 0(4^3.) 

2(aL + aR) 2ala^p,[aL + aR) [ai + a/?) 

Using Laplace's method [56] to asymptotically evaluate ( IHOj) . the substitution of (!82|) and (!83|l 
into ([79]) yields 



sgn(2;(s)) sgn(2;(u))) du ~'~ ^^■'^ 



JO ala%{aL + aR) 

- ( ^ \ + 2 2 / ^ ^ \ aHPe(0,t;|0) + -^(0,^|0) U^ 

By applying Lemma H] and simplifying we finally arrive at ( l56l) . □ 



19 



References 



[1] M. di Bernardo, C.J. Budd, A.R. Champneys, and P. Kowalczyk. Piecewise- smooth Dynam- 
ical Systems. Theory and Applications. Springer- Verlag, New York, 2008. 

[2] A.J. Van der Schaft and J.M. Schumacher. An Introduction to Hybrid Dynamical Systems. 
Springer- Verlag, New York, 2000. 

[3] R.I. Leine and H. Nijmeijer. Dynamics and Bifurcations of Non-smooth Mechanical Systems., 
vohimc 18 of Lecture Notes in Applied and Computational Mathematics. Springer- Verlag, 
Berlin, 2004. 

[4] S. Banerjcc and G.C. Verghese, editors. Nonlinear Phenomena in Power Electronics. IEEE 
Press, New York, 2001. 

[5] Z.T. Zhusubaliyev and E. Mosekilde. Bifurcations and Chaos in Piecewise-Smooth Dynam- 
ical Systems. World Scientific, Singapore, 2003. 

[6] T. Puu and I. Sushko, editors. Business Cycle Dynamics: Models and Tools. Springer- Verlag, 
New York, 2006. 

[7] H. Dankowicz and A.B. Nordmark. On the origin and bifurcations of stick-slip oscillations. 
Phys. D, 136:280-302, 2000. 

[8] M.H. Frcdriksson and A.B. Nordmark. On normal form calculation in impact oscillators. 
Proa. R. Soc. A, 456:315-329, 2000. 

[9] M. di Bernardo, C.J. Budd, and A.R. Champneys. Normal form maps for grazing bifurcations 
in n-dimensional piecewise-smooth dynamical systems. Phys. D, 160:222-254, 2001. 

[10] N. Berglund and B. Gcntz. Noise-Induced Phenomena in Slow-Fast Dynamical Systems. 
Springer, New York, 2006. 

[11] B. Lindner, J. Garcia-Ojalvo, A. Neiman, and L. Schimansky-Geier. Effects of noise in 
excitable systems. Phys. RepoHs, 392:321-424, 2004. 

[12] A.S. Pikovsky and J. Kurths. Coherence resonance in a noise-driven excitable system. Phys. 
Rev. Lett., 78(5):775-778, 1997. 

[13] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni. Stochastic resonance. Rev. Modern 
Phys., 70(l):223-287, 1998. 

[14] T.C.L. Griffin. Dynamics of Stochastic Nonsmooth Systems. PhD thesis. University of 
Bristol, 2005. 

[15] R. Wackerbauer. Noise-induced stabilization of one-dimensional discontinuous maps. Phys. 
Rev. E, 58(3):3036-3044, 1998. 



20 



[16] L. Zhang, P. Shi, C. Wang, and H. Gao. Robust filtering for switched hnear discrete-time 
systems with polytopic uncertainties. Int. J. Adapt. Control Signal Process., 20:291-304, 
2006. 

[17] W. Zhang, J. Hu, and J. Lian. Quadratic optimal control of switched hnear stochastic 
systems. Syst. Contr. Lett., 59:736-744, 2010. 

[18] P.H.E. Tiesinga. Precision and rehabihty of periodically and quasiperiodically driven 
integrate-and-fire neurons. Phys. Rev. E, 65(4):041913, 2002. 

[19] T. Griffin and S. Hogan. Dynamics of discontinuous systems with imperfections and noise. 
In G. Rega and F. Vestroni, editors, lUTAM Symposium on Chaotic Dynamics and Control 
of Systems and Processes in Mechanics., pages 275-285. Springer, 2005. 

[20] M.F. Dimentberg and D.V. lourtchenko. Random vibrations with impacts: A review. Non- 
linear Dyn., 36:229-254, 2004. 

[21] M.F. Dimentberg and A.I. Menyailov. Response of a single-mass vibroimpact system to 
white-noise random excitation. Z. Angew. Math. Mech., 59(12):709-716, 1979. 

[22] M. Fogli, P. Bressolette, and P. Bernard. The dynamics of a stochastic oscillator with 
impacts. Eur. J. Mech. A-Solids, 15(2):213-241, 1996. 

[23] N Sri Namachchivaya and J.H. Park. Stochastic dynamics of impact oscillators. J. Appl. 
Mech. Trans. ASME, 72(6):862-870, 2005. 

[24] J. Feng, W. Xu, H. Rong, and R. Wang. Stochastic responses of Duffing- Van der Pol vibro- 
impact system under additive and multiplicative random excitations. Int. J. Non-Linear 
Mech., 44:51-57, 2009. 

[25] P.D. Christofides and N.H. El-Farra. Control of Nonlinear and Hybrid Process Systems. 
Designs for Uncertainty, Constraints and Time-Delays. Springer, New York, 2005. 

[26] W. Feng and J.-F. Zhang. Stability analysis and stabilization control of multi-variable 
switched stochastic systems. Automatica, 42:169-176, 2006. 

[27] P. Mhaskar, N.H. El-Farra, and P.D. Christofides. Robust hybrid predictive control of 
nonlinear systems. Automatica, 41:209-217, 2005. 

[28] B. Lindner and L. Schimansky-Geier. Coherence and stochastic resonance in a two-state 
system. Phys. Rev. E, 61(6):6103-6110, 2000. 

[29] D.J.W. Simpson and R. Kuske. Mixed-mode oscillations in a stochastic piecewise-linear 
system. Phys. D, 240:1189-1198, 2011. 

[30] A.F. Filippov. Differential Equations with Discontinuous Righthand Sides. Kluwer Academic 
Publishers., Norwell, 1988. 

[31] M. Wiercigroch and B. De Kraker, editors. Applied Nonlinear Dynamics and Chaos of 
Mechanical Systems with Discontinuities., Singapore, 2000. World Scientific. 



21 



[32] B. Brogliato. Nonsmooth Mechanics: Models, Dynamics and Control. Springer- Ver lag, New 
York, 1999. 

[33] B. Blazejczyk-Okolewska, K. Czolczynski, T. Kapitaniak, and J. Wojewoda. Chaotic Me- 
chanics in Systems with Impacts and Friction. World Scientific, Singapore, 1999. 

[34] J. Awrejcewicz and C. Lamarque. Bifurcation and Chaos in Nonsmooth Mechanical Systems. 
World Scientific, Singapore, 2003. 

[35] R.A. Ibrahim. Vibro-Impact Dynamics., volume 43 of Lecture Notes in Applied and Com- 
putational Mechanics. Springer, New York, 2009. 

[36] C.K. Tse. Complex Behavior of Switching Power Converters. CRC Press, Boca Raton, FL, 
2003. 

[37] A.F. Filippov. Differential equations with discontinuous right-hand side. Mat. Sb., 51(93):99- 
128, 1960. English transl. Amer. Math. Soc. Transl. 42(2):199-231, 1964. 

[38] P. Casini, O. Giannini, and F. Vestroni. Experimental evidence of non-standard bifurcations 
in non-smooth oscillator dynamics. Nonlinear Dyn., 46(3):259-272, 2006. 

[39] A.C.J. Luo and B.C. Gegg. Stick and non-stick periodic motions in periodically forced 
oscillators with dry friction. J. Sound Vib., 291:132-168, 2006. 

[40] M. Johansson. Piecewise Linear Control Systems., volume 284 of Lecture Notes in Control 
and Information Sciences. Springer- Verlag, New York, 2003. 

[41] M. di Bernardo, K.H. Johansson, and F. Vasca. Self-oscillations and sliding in relay feedback 
systems: Symmetry and bifurcations. Int J. Bifurcation Chaos, 11(4):1121-1140, 2001. 

[42] Z. Schuss. Theory and Applications of Stochastic Differential Equations. Wiley, New York, 
1980. 

[43] C.W. Gardiner. Handbook of Stochastic Methods for Physics, Chemistry and the Natural 
Sciences. Springer- Ver lag. New York, 1985. 

[44] J. Grasman and O.A. van Herwaarden. Asymptotic Methods for the Fokker-Planck Equation 
and the Exit Problem in Applications. Springer, New York, 1999. 

[45] Ya.Z. Tsypkin. Relay Control Systems. Cambridge University Press, New York, 1984. 

[46] G.F. Franklin, J.D. Powell, and A. Emami-Naeini. Feedback Control of Dynamic Systems. 
Prentice Hall, Upper Saddle River, NJ, 2002. 

[47] R.C. Dorf and R.H. Bishop. Modern Control Systems. Prentice Hall, Upper Saddle River, 
NJ, 2001. 

[48] K.J. Astrom and R.M. Murray. Feedback Systems. An Introduction for Scientists and Engi- 
neers. Princeton University Press, Princeton, NJ, 2008. 



22 



[49] K.H. Johansson, A. Rantzer, and K.J. Astrom. Fast switches in relay feedback systems. 
Automatica, 35:539-552, 1999. 

[50] K.H. Johansson, A.E. Barabanov, and K.J. Astrom. Limit cycles with chattering in relay 
feedback systems. IEEE Trans. Automat. Contr., 47(9):1414-1423, 2002. 

[51] Y. Zhao, J. Feng, and C.K. Tse. Discrete-time modeling and stability analysis of periodic 
orbits with sliding for switched linear systems. IEEE Trans. Circuits Systems I Fund. Theory 
Appl, 57(ll):2948-2955, 2010. 

[52] M. Tanelli, G. Osorio, M. di Bernardo, S.M. Savaresi, and A. Astolfi. Existence, stability 
and robustness analysis of limit cycles in hybrid anti-lock braking systems. Int. J. Contr., 
82(4):659-678, 2009. 

[53] M. di Bernardo, K.H. Johansson, U. Jonsson, and F. Vasca. On the robustness of periodic 
solutions in relay feedback systems. In 15th Triennial World Congress, Barcelona, Spain, 
2002. 

[54] H. Risken. The Fokker-Planck Equation: Methods of Solution and Applications. Springer- 
Verlag, New York, 1984. 

[55] Z. Schuss. Theory and Applications of Stochastic Processes. Springer, New York, 2010. 

[56] CM. Bender and S.A. Orszag. Advanced Mathematical Methods for Scientists and Engineers. 
International Series in Pure and Applied Mathematics. McGraw-Hill, New York, 1978. 

[57] I. Karatzas and S.E. Shreve. Trivariate density of Brownian motion, its local and occupation 
times, with apphcation to stochastic control. Ann. Prob., 12(3):819-828, 1984. 

[58] I.V. Girsanov. On transforming a certain class of stochastic processes by absolutely contin- 
uous substitution of measures. Theory Prob. Appl, 5(3):285-301, 1960. 

[59] B. Oksendal. Stochastic Differential Equations: An Introduction with Applications. Springer, 
New York, 2003. 

[60] 1. Karatzas and S.E. Shreve. Brownian Motion and Stochastic Calculus. Springer, New 
York, 1991. 

[61] V.E. Benes, L.A. Shepp, and H.S. Witsenhausen. Some solvable stochastic control problems. 
Stochastics, 4:39-83, 1980. 

[62] Z. Qian and W. Zheng. Sharp bounds for transition probability densities of a class of 
diffusions. CR. Acad. Sci. Paris, Ser. I, 335:953-957, 2002. 

[63] Z. Qian, F. Russo, and W. Zheng. Comparison theorem and estimates for transition proba- 
bility densities of diffusion processes. Probab. Theory Relat. Fields., 127:388-406, 2003. 

[64] M. Gradinaru, S. Herrmann, and B. Roynette. A singular large deviations phenomenon. 
Ann. I. H. Poincare, 37(5):555-580, 2001. 



23 



[65] W. Zhang. Transition density of one-dimensional diffusion with discontinuous drift. IEEE 
Trans. Automat. Contr., 35(8):980-985, 1990. 

[66] A. Colombo, M. di Bernardo, S.J. Hogan, and M.R. Jeffrey. Bifurcations of piecewise smooth 
flows: perspectives, methodologies and open problems. Submitted to: Phys. D, 2012. 



24 



